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The perturbation of a light field might afi^ect preheating and hence generate a contribution to the 
spectrum and non-gaussianity of the curvature perturbation i^. The field might appear directly in the 
preheating model (curvaton-type preheating) or indirectly through its effect on a mass or coupling 
(modulated preheating). We give general expressions for C, based on the 5N formula, and apply 
them to the cases of quadratic and quartic chaotic infiation. For the quadratic case, curvaton-type 
preheating is ineffective in contributing to but modulated preheating can be effective. For quartic 
inflation, curvaton-type preheating may be effective but the usual SN formalism has to be modified. 
_ ^ We see under what circumstances the recent numerical simulation of Bond et al. [0903.3407] may 

^ , be enough to provide a rough estimate for this caseQ 

PACS numbers: 98.80.Cq 

I. INTRODUCTION 

As far as we can tell, the observed cosmological perturbations originate from a primordial 
curvature perturbation C(x), that is present at the epoch when the shortest cosmological 
scale approaches horizon entry. That epoch corresponds to a temperature of order 1 MeV 
and an age t ^ la. At this stage the curvature perturbation is time-independent because 
CIh' the cosmic fluid is radiation-dominated to very high accuracy, but it may be time dependent 

^ . at earlier times. 

It is thought that C(^:^) each position is determined by the values (?!>*(x), of one or 
more scalar fields evaluated at some initial epoch during inflation. (We shall consistently 
'sj" ■ use a star to denote quantities evaluated at this epoch.) The corresponding perturbations 

^ \ 5(j)* are supposed to be generated from the vacuum fluctuation, which requires the fields to 

be light. 

According to the original paradigm, ^ is generated entirely from the perturbation Scj):, of 
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. the inflaton field in a slow roll inflation model. In that case, C is present already at the 

initial epoch, being constant thereafter. The non-gaussianity of C is in this case very small, 
and almost certainly undetectable. 

According to alternative paradigms, the perturbation Sx* of a light field other than the 
I inflaton generates a significant (maybe dominant) contribution. Such a contribution is ini- 

tially negligible, growing to its final value later, and in general only after inflation is over. 
The non-gaussianity in this case can be detectable. 

In this paper we consider the growth of C that may occur during preheating. Preheating 
' may be defined as the loss of energy by time-dependent scalar fields, through mechanisms 

5^ , other than single particle decay. Preheating is typically followed by reheating, with an 

intervening era between the two. Reheating may be defined as the practically complete 
thermalization of a gaseous cosmic fluid, including the case that the constituents of the gas 
correspond to an oscillating scalar field. 

The simplest mechanism for preheating is the potential 

F=imV + ^5Vx^ (1) 

with (j) the inflaton field, and x = during inflation so that we deal with the quadratic 
('chaotic') inflation potential. If no other terms of the action are relevant, (f> oscillates and 
loses energy by creating x particles through what is called parametric resonance 

UM- As 

in this example, the field x is usually different from the field 0, but parametric resonance 
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works in just the same way if they are the same. It also works with different forms of the 
potential, such as Q 

V^\m' + \9W, (2) 

with again the inflaton field. (That case is often called massless preheating.) Given a po- 
tential that allows preheating, it can happen that the produced field x can decay sufficiently 
rapidly into (say) a pair of fermions, so that 4> can lose its energy before even one oscillation 
takes place. That is called instant preheating Q. There is also what is known as tachyonic 
preheating [8|, which in its simplest form invokes the potential 

V{x) ^Vo- \m\^ + -^\x\ (3) 

Tachyonic preheating occurs because the effective mass-squared —w? is initially to be posi- 
tive. During the sign change of rr? the vacuum fluctuation of x to be promoted to a classical 
perturbation which is then amplified as x rolls off the hilltop. Parametric resonance can 
take place is x passes through the minimum, followed by more amplification of the original 
perturbation and so on. 

These basic preheating mechanisms have been considered within several different scenar- 
ios. In most of them, preheating takes place immediately after inflation so that the oscillating 
field is the one involved in the inflation model. (If the inflation model is hybrid or multi-field 
there are two or more of these fields.) In some scenarios though, preheating takes place after 
a later phase transition and the oscillating field played no role during inflation. 

As we will explain, any growth of Q occurring during preheating may persist for some time 
afterwards. The growth will generally terminate at or before reheating. (If there is further 
growth during or after reheating then that growth should be treated as a separate process.) 
For clarity we will just talk about 'growth during preheating' on the understanding that any 
subsequent growth prior to reheating is to be included in the discussion. 

Regarding its ability to generate a contribution to C, preheating has some points of simi- 
larity with reheating. Two distinct mechanisms exist in the reheating case. In the original 
'curvaton' mechanism*-'^ (Tol - [l3 |. the light field responsible for generating the contribution 
to C is the oscillating scalar field, whose decay is responsible for reheating. Later, the 'mod- 
ulated reheating' mechanism was proposed There, the relevant light field acts only 
indirectly, by affecting the decay rate of the inflaton so that we have r(x(x)) where that F 
is the decay rate and x is the perturbed light fleld. With this in mind, we can distinguish 
two versions of the preheating scenario. In a 'curvaton-type preheating' scenario the light 
field is directly involved in preheating; either a field created by preheating, or an oscillating 
field that is responsible for the creation. In a 'modulated preheating' scenario the light field 
instead acts indirectly, by affecting a preheating parameter such as the coupling g in Eq. P]). 
In that case we have g{x) where x is the light field. 

Modulated preheating or reheating can be implemented within any scenario, but becomes 
predictive only when one specifies the dependence of the parameter on the field (ie. the func- 
tion F(x) or gijif). Curvaton reheating and curvaton-type preheating are more restricted, 
but more predictive because F, g and so on are taken to be constants. We shall show that 
curvaton-type preheating cannot occur with the potential ([T]) if it is supposed to hold also 
during infiation. 

What we are calling modulated preheating has not so far been studied, though its possi- 
bility has been recognized [iT]]. (For related works see [3.) The possibility of what we are 
calling curvaton-type preheating has been considered in a few papers. The papers [l9ll20j. 
consider the quadratic inflaton potential ([IJ, taking the potential to be valid also during 
inflation. The papers j2()hi22i] consider tachyonic preheating. The papers l,23J c onsider a 



multi-field generalization of Eq. ([T]) with instant preheating. The papers [24| - |27l | consider 



The mechanism was also proposed in the context of a bouncing universe; see 0| for a discussion and earUer 
references. 
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the massless preheating scenario of Eq. ([2]) . Reference [2g] considers preheating at the end 
of a hybrid inflation model involving the Hig gs field. In these works, preheating starts as 
soon as inflation is over. In another work [23 the curvaton field causes preheating, when it 
starts to oscillate long after inflation is over. In [s^ preheating is caused by the oscillation 
of a flat direction of the MSSM, again long after inflation is over. 

Most of the papers cited above use cosmological perturbation theory, at either first order 
mil mill or second order [H-ill,!!!. We instead use the 6N formalism fallH,!!!. We 
begin in Section [II] by recalling the basic properties of the curvature perturbation. In Section 
lllll we carefully set up the SN formalism. In Section HVl we show that the quadratic potential 
Eq. ([IJ cannot provide curvaton-type preheating, contrary to what was assumed in [Tgl. [20j. 
In Section |V] we consider modulated preheating with the quadratic potential. In Section 
IVII we consider curvaton-type preheating the quartic potential ([2]) (massless preheating). In 
contrast with other applications of the 6N formalism, the perturbations of the light fields 
generated after cosmological scales leave the horizon are likely to be significant in this case. 
Ignoring them, we show how to use the recent numerical simulation of [27|| to give an order of 
magnitude estimate of /nl with the unperturbed value of x* set equal to zero. It disagrees 
with the estimate obtained in [25] using second order cosmological perturbation theory. In 
a concluding section we summarise our finding, and point to future directions for research. 
In two appendices we extend the discussion of Section |IVl 



II. THE CURVATURE PERTURBATION 



A. Definition and evolution 



The non-perturbative definition of the primordial curvature perturbation C is described 



in for instance 37[, where original references can be found. The components of the metric 
tensor are smoothed on a comoving scale R and one considers the super-horizon regime 
aR > where H = a/a is the Hubble parameter and a{t) is the scale factor normalised 
to 1 at present. "^^^ The energy density p and pressure P are smoothed on the same scale. 
One considers the slicing of spacetime with uniform energy density. The spatial metric is 
written as 

g.,(x,i)^a2We^^(-'*)(/enx,0),^., (4) 

where h is traceless so that has unit determinant. The smoothing scale is chosen to 
be somewhat shorter than the scales of interest, so that the Fourier components of C, on 
those scales is unaffected by the smoothing. The threading of spacetime is taken to be 
orthogonal to the slicing. The time dependence of the locally defined scale factor a(x, = 
a(i)exp(^) defines the rate at which an infinitesimal comoving volume V expands: V/V = 
3d(x, t)/a(x, t). 

Under the reasonable assumption that the Hubble scale is the biggest relevant dis- 
tance scale, the energy continuity equation diVp) = —PdV at each location is the same as 
in a homogeneous universe; as far as the evolution of p is concerned, we are dealing with 
a family of separate homogeneous universes. With the additional assumption that the ini- 
tial condition is set by scalar fields during inflation (adopted here), the smoothed hij{t) is 
time-independent after smoothing and then the separate universes are homogeneous as well 
as isotropic. 

Since we are working on slices of uniform p, the energy continuity equation can be written 
p{t) = -3 \H{t) + C(x, t)] [p{t) + P(x, t)] . (5) 



Smoothing a function g{x) means that g at each location is replaced by its average within a sphere of 
coordinate radius R around that position. The averaging may be done with a smooth window function 
such as a gaussian. The smoothed function is supposed to have no significant Fourier components with 
coordinate wavenumbers k ^ R, which means that its gradient at a typical location is at most of order 
1/R. A function g with that property is said to be 'smooth on the scale R\ 
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One write 

P{x,t)^P{t)+dP{x,t), (6) 

so that 6P is the pressure perturbation on the uniform density shces, and choose P{t) so 
that the unperturbed quantities satisfy the unperturbed equation p = —3H{p + P). Then 

C = (7) 

^ P + P + 6P ^' 

This gives C if we know p{t) and P(x, t). It makes C time-independent during any era when 
P is a unique function of p [s^ , Issj (hence uniform on shces of uniform p) . The pressure 
perturbation is said to be adiabatic in this case, otherwise it is said to be non-adiabatic. 

The key assumption in the above discussion is that in the superhorizon regime certain 
smoothed quantities (in this case p and P) evolve at each location as they would in an 
unperturbed universe. In other words, the evolution of the perturbed universe is that of a 
family of unperturbed universes. This is the separate universe assumption, that is useful 
also in other situations [37| . 

The primordial curvature perturbation ^ is directly probed by observation on 'cosmological 
scales' corresponding to roughly e~^^HQ^ < < H^^. These scale begin to enter the 
horizon when when T ^ 1 MeV. The Universe at that stage is radiation dominated to very 
high accuracy, implying P = p/3 and a constant curvature perturbation which we denote 
simply by C(^)- When cosmological scales are the only ones of interest, one should choose 
the smoothing scale as i? ~ e'^^H^^ . Unless stated otherwise, we make this choice. 

Within a given scenario, C will exist also on smaller inverse wavenumbers, down to some 
'coherence length' which might be as low as ^ e~^^HQ^ (the scale leaving the horizon 
at the end of inflation). If one is interested in such scales, the smoothing scale R should be 
chosen to be (somewhat less than) the coherence length. 



B. Correlators 



Theories will predict the correlators of cosmological perturbations. The two-point corre- 
lator of a perturbation (/(x) is written {g{xi)g{x2)) and so on for higher correlators. The 
bracket basically denotes an ensemble average, with our universe a typical member of the 
ensemble. If the perturbations originate as a vacuum fluctuation during inflation, the en- 
semble average is the vacuum expectation value. If the universe is homogeneous this makes 
the correlators translation invariant, which we will assume. 

Given translation invariance, the ergodic theorem holds whereby the bracket can be taken 
to be a spatial average. Thus {g{xi)g{x2)) is the spatial average of g{xi)g(ii2) at fixed 
xi — X2, within any volume with size L much bigger than |xi — X2I. For any finite L, 
the volume average differs from the the ensemble average. The difference is called cosmic 
variance. 

To exploit the translation invariance it is convenient to work with Fourier components g^. 
Using a box of coordinate size L we have 

j d^xg{^)e-'^ ''. (8) 

L-^ J2 ffke'" " - (2^)"' / d'kgy^e'''-^. (9) 

Since the box imposes a periodic boundary condition, only the wavenumbers kL ^ 1 are 
physically significant. For them one can use the final equality. According to the ergodic 
theorem, one can regard (gkiffk2) SI'S an average over the points within a cell of momentum 
space, up to the cosmic variance. 



5k = 
5(x) = 



The quantity (g{x)) is then just a number which can be set to zero by choice of the unperturbed quantity. 
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The dependence of (the loop contributions to) the correlators upon the box size, reflects 
the fact that the correlators are in position space are spatial averages within the box. Ob- 
servations are available within the observable universe, with coordinate size of order Except 
for the low multipoles of the CMB, all observations probe scales k ^ Hq. To handle them, 
one should choose the box size as L = ^^o"^ [22]- A smaller choice would throw away some 
of the data while a bigger choice would make the spatial averages unobservable. 

Low multipoles of the CMB anisotropy explore scales of order H^^ /i not very much 
smaller than H^"^ . To handle them one has to take L bigger than H^^ . For most purposes, 
one should use a box, such that ln(Li?o) is just a few (ie. not exponentially large) [33, Ho, |4l| . 

To summarise, the box size should normally be taken as L = H^^. The exception is the 
case of low CMB multipoles, where one should normally take L ^ H,^^ with In(iiJo) not 
exponentially large. A box satisfying these requirements will be called a minimal box. 

The spectrum Pg is defined in terms of the two-point correlator by 

(5k5k') = (27r)=^<53(k + k')P<,(fc). (10) 
One also defines Vg{k) = (fc^/27r^)Pg(fc), also called the spectrum. An equivalent expression 

m is 

(5(y).9(x + y)) = (27r)-3 J d?kPg{k)e^'^--, (11) 

where the bracket can be regarded as the spatial average with respect to y. 

The mean-square perturbation, evaluated within a box of size L and smoothed on scale 
R, is 

f^"' dk 

{g\^))=j^_^ -Vg{k)^VgHL/R). (12) 

The final estimate is valid if Pg{k) is more or less scale- independent as will usually be the 
case for primordial perturbations. The spectrum of dg/dxi (the derivative with respect to 
one of the coordinates) is kfVg, hence the mean-square gradient is 

((|Vff|)2) = J^_^ fk^Vgik) ^Vg/R' ^ {g')/R\ (13) 

For a gaussian perturbation, the spectrum defines all correlators. The n-point correlators 
vanish for odd n, while for even n 

(.9ki5k2.9k3.gk4> = (.9ki.9k2)(5k3 5k4> + (.9ki.9k3)(5k3 5k2) + (gkigk4)(5k2 5k3), (14) 

and so on. For a non-gaussian contribution one has to specify more quantities, starting with 
the bispectrum Bg: 

(5kiS'k2<7k3) - (277)3^3 (ki -I- k2 -)-k3)Sg(/ci,fc2,fc3) (15) 

Bg{ki,k2,ks,) = Sg(fci, ^2, fcs) [Pg(A:i)Pg(fc2) -f cyclic permutations] . (16) 

The quantity Bg is called the reduced bispectrum. For C, one uses /nl = (5/6)S^. 

From observation of the CMB anisotropy and the galaxy distribution, we know that V^ik) 
is almost scale-independent with the value (5 x 10~^)^ ~ 10~^. Also, we know that C, is 
gaussian to high accuracy. Taking /nl to be practically scale independent ('local' form) as 
is predicted by the models that we consider, the bound at 95% confidence level according 
to the first of [H is -4 < /nl < 80. 

III. THE SN FORMULA 
A. The general formula 



Consider now, the analogue of Eq. (U) for a generic slicing: 



(17) 
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where 7 has unit determinant. Take the initial shce to be flat (meaning; that -0 = 0) and the 
final slice to be of uniform energy density so that = C,. This gives [32l | 

C(x,t) = <57V(x,t), (18) 

where A'' is the number of e-folds of expansion between the initial and final slice. It is 
independent of the choice of the initial flat slice because the expansion going from one flat 
slice to another is uniform. By virtue of the separate universe assumption, this formula 
allows one to calculate C(x, t) given the evolution of the scale factor in some family of 
unperturbed universes. 

To proceed, we invoke light fields, taken to be canonically normalized. Ignoring mixed 
derivatives for simplicity, a light field is defined [s^] as one whose effective mass-squared 
mf during inflation satisfles 

|m2|«i/^ m^ = d^V/d(l)l (19) 

As each scale k leaves the horizon during inflation, the vacuum fluctuation of each light 
field is converted to a classical perturbation [s^, 1131 • At a given epoch during infiation, the 
classical fields 0i(x) are therefore smooth on the horizon scale. 

Now we take the initial epoch to be the one when the smoothing scale R leaves the horizon, 
and make a crucial assumption. At each location, one or more of the light fields provides the 
initial condition for the evolution of the local scale factor (along with unperturbed quantities 
including parameters of the field theory and the values of any relevant unperturbed fields). 
Then we have [H] 

C(x, t)=SN = N{p{t), 0t (x), (/.^(x), • • •) - N{p{t), 0t , • • •), (20) 

where the star denotes the field values at the initial epoch. Since the typical magnitude of 
a light field perturbation is of order i/*, we need values of ^^(x) in a range A(/)* ^ (i7*/27r), 
centred on the unperturbed values . 

These initial values (/)*(x) are smooth on the scale R, and so are the initial values of 
p and P. If the subseqent evolution of the light fields were classical, they would remain 
smooth on the scale i?, and so would p and P. In fact, the vacuum fluctuation continues 
to generate classical fleld perturbations, so that by the end of inflation they are present 
on scales iaH)~^^ < < R. In all known cases except that of massless preheating, 
these smaller scale perturbations in the light field have little effect on the perturbations 
in p and P on scales of interest k^^ > R and so can be ignored. Then, by virtue of the 
separate universe assumptions, the evolution of p and P at each location is the same as 
in an unperturbed universe, and so is the evolution of N. We can evaluate the curvature 
perturbation by considering a family of unperturbed universes! This is the usual version of 
the SN formalism, which we use except when dealing with massless preheating. 



B. Power series in the field perturbations 

For each field we can write 

</>*(x) = (/.:+<50*(x). (21) 

In all cases so far considered, we can expand N {p{t) , (pKx) , (p2{^) , ■ ' ') &s a low-order power 
series about the unperturbed field values: 

c(x, t) = j2 mm w + ^ E ^'j- (^)^<^* (^)'^<^.* (x) + . . . , (22) 

where a subscript i denotes d/dcj)* evaluated at the unperturbed point of field space. To get 
the observed quantity C(x) this should be evaluated after the Ni{t) and so on settle down 
to their final values, which we denote simply by Ni and so on. 

This formula was first given in [s^l following the non-linear proof in [s^l of ^ = 5N . To 
first order in the fleld perturbations, the formula and the proof were given in [3l| . which 
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focuses on the case that C is generated during inflation. Related formulas for the evolution 



of during inflation were given in [431] (first order) and in 4J| (non-linear). 

The light field perturbations have completely negligible non-gaussianity [46|, and their 
spectrum soon after horizon exit for the scale k is given by [37j 

Vsc^Ak)-{Hk/2^f, (23) 

where Hk can be taken as the value of H at horizon exit. The observed near-gaussianity 
of C requires that Eq. (|22p is well approximated by one or more linear terms, with the field 
perturbations nearly gaussian. Taking the initial epoch to be just after horizon exit we then 
have 



n(fc)-E^'(^) ■ (24) 

Assuming that the reduced bispectrum /nl is big enough to be observable (say |/nl| ^1), 
it is generated by the quadratic terms, the slight non-gaussianity of the 5(j)i being negligible. 
In this case, /nl is almost scale-independent (local form). 



C. The inflaton contribution to C, 



We will suppose that there is standard slow-roll inflation, which is effectively single-fleld 
because only the inflaton field cj) has significant variation during inflation. 

The slow-roll approximation gives V ~ 3MpH^ and ~ —V'/3H, where V{(f>) is the 
potential during slow-roll inflation. This gives the well-known relation 

N{k) = A/p 2 / (25) 

where N{k) denotes the number of e-folds of slow-roll inflation occurring after the epoch of 
horizon exit for the scale k~^. The inflationary potential will determine the field value </>end 
at which slow-roll infiation ends. 

One usually denotes N{Ho) simply by N. It is determined by the inflationary energy scale 
i?* and the post-inflationary evolution of the scale factor. A typical value, which we adopt, 
is iV ~ 55. For a smaller scale, N{k) = N — \n{k/Ho). For the smallest cosmological scale 
(enclosing mass lO^M© or so) N{k) ~ N — 15 = 40. The initial epoch, when the value 
of (p is denoted by 0*, can be soon after that. The important feature of Eq. (I^5|) is that 
it determines the unperturbed field value and hence Vq^ , given the inflationary potential 
and the post-inflationary cosmology. 

Perturbing </) with any other light fields fixed gives a constant contribution C0(x) to C: 
because it represents a shift back and forth along the trajectory which does not alter the 
subsequent relation between p and P. The perturbation (50* has typical magnitude if,/27r ^ 
(/)*. Since C,^ represents a shift along the inflaton trajectory, A^^ — —H/cj) evaluated at the 
initial epoch. Using the slow-roll approximation one finds that is practically linear, 
making it practically gaussian j45| . Taking the initial epoch soon after horizon exit, 

Vc,{k) =Nl{Hk /27rf. (26) 

If 7-"^^ dominates, the spectral tilt is n — 1 = 2?7 — 6e, where 77 = MpV" /V and e = 
M^{V' /V)^ /2 evaluated at horizon exit. 



D. The non-inflaton contribution 



Now we consider the contribution of light fields other than the inflaton. We suppose that 
only one of them contributes significantly, and call it a. Its contribution ^^{'x.,t) is initially 
negligible, but it can grow by generating a non-adiabatic pressure perturbation. This might 
happen during any era, except one of complete matter domination (P = 0) or complete 
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radiation domination (P = p/3). Possibilities in chronological order include generation 
(i) during multi- field inflation (ii) at the end of inflation [13], (iii) during preheating 
as we discussed in the Introduction, (iv) during modulated reheating (la - flTl l48l | , (v) at a 
modulated phase transition ^49*] and (vi) before a second reheating through the curvaton 
mechanism |14i] - In this paper we focus on preheating. 

Wc interested in the value of Co- after it settles down to its final constant value, which we 
denote simply by C(t(x) 

C.(x) ^ C(x)-C0(x) (27) 
= ^A^(a,(x),0,) + (5[iV(a,(x),0,(x)))-7V(cr,(x),0,)-7V(CT,,0,(x))]. (28) 

We make the foUowign approximation, which is usually adequate: 

C,(x) ~ 5iV(a4x),0,) (29) 

~ N^Sa^x) + ^N^^ {Sa4x)f . (30) 

Eq. ([^ ignores mixed derivatives d"^ N / d(j)^,d<7^ and so on, while Eq. (15(11) ignores d^N/da^ 
and higher derivatives. 

The unperturbed value cr* is defined as the spatial average of ct(x) within the minimal box. 
It is in general a free parameter whose value has to be specified along with the parameters 
of the Lagrangian. To get some idea about is likely value though, we can use the fact that 
the mean square perturbation, smoothed on scale Li and evaluated within a larger box of 
size L2, is 

{5a,Y^ 1 ^n..(fc)^ (i?*/27r)2ln(L2/ii), (31) 

where in the final equality we set Va, ^ {H^,/2tt)'^. To use this formula we can take Li = L 
with L a minimal box centred on the observable universe. Then the formula says that the 
unperturbed value of cr* (ie. it's spatial average within the minimal box) will be at least 
of order i/*/27r if the observable universe is located at a typical position. If there were a 
very large number of e-folds of almost-exponential inflation before the observable universe 
left the horizon, we can go further; taking the formula literally we can choose lii{L2 / Li) to 
be much bigger than 1, and conclude that the value of the unperturbed fleld for a typical 
location of our universe will be much bigger than H^,/2n. (As ln(L2/£i ) is increased the 
calculation leading to Eq. will at some stage become out of control [40| but we don't 
need an enormous value of ln(L2/^i) to arrive at the above conclusion.) Of course, it might 
be that we live in a special place where the the unperturbed x* is much less than H^/2tt. 

Using Eqs. dH), (Uni), (HH) and and taking Sa^ to be gaussian with a scale- independent 
spectrum Vsa, ~ (^*/27r)^, we can evaluate the spectrum and bispectrum of the curvature 
perturbation. One finds 

Vc^ = + 75'°°? (32) 

ptrec ^ ^2p^^^ ^ NUH,/2Trf (33) 

pio;P(A:) = ^-NlV^s.Mk) = NlHkL)(^^y , (34) 



/NL 



/i^r+zirr (35) 
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f^T = ^ ( ) ^ (36) 



and 



6 V n / 

f°l^\kuk2,ks) ^ ln{kL)^(^^j . (37) 

The labels 'tree' and 'loop' refer to the Feynman-like diagrams of ^Ed\. The loop contribu- 
tions correspond to formally divergent integrals, but they are regularized by setting Vq^ = 
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for k ^ L~^, with L the box size. Eq. (p4)) was given in [si*]. It is vahd for k ^ L~^, which 
is anyhow required so that one can ignore the artificial periodicity of the Fourier series. 
Eq. (|36p was first given in [1]. Eq. (|37p was given in [52l |. with fc^ taken to have a common 
value k 3> to rough order of magnitude. 
We see that 

= ln(fcL)(|l^)Vr' (38) 

- ln(A:L) (|7^)'/Nr. (39) 

By examining the integral 52] leading to Eq. (P7)) , one can see that Eq. (15^)1 becomes exact 
in the squeezed configuration k = ki <^ k2 — k^. This has not been noticed before, and is 
quite interesting because the observational constraint on a local /nl comes mostly from the 
squeezed configuration [42| . 

We see that the tree contributions dominate if 

while the loop contributions dominate in the opposite case. An equivalent statement is that 
the tree contributions dominate if the linear term of Eq. pop dominates, while the loop 
contributions dominate if the quadratic term of Eq. (pO|) dominates. (The statements are 
equivalent because the mean-square of (5ct*, smoothed on the scale fc, is ln(/cL)(i?*/27r)^.) 

When comparing the loop contribution with observation one should normally set L ~ Hq^ , 
except for the low CMB multipoles where one should choose L ^ Hf^^ with In(fci) ^ 1. 
With the choice L = Hq^, In(fcL) ~ 5 for the scales explored by the CMB multipoles with 
£ ~ 100, while In(fcL) ~ 10 for the scales explored by galaxy surveys. This increase would 
be observable [39]. 

If Co- ~ C dominates, the near-gaussianity of C requires [1^ that the tree contributions 
dominate. Then, by taking the initial epoch to be soon after horizon exit, one finds spectral 
tilt 

n - 1 = 2?7^^ - 2e, (41) 

where t?^^ = Ml.{d^V/dal)/V 

E. Evolution of a 

As noted earlier, we can usually take the relevant light fields 4>i{'x., t) to be smooth on the 
scale R even after the initial epoch. Then we can invoke the separate universe assumption, 
making their evolution at each location the same as in some unperturbed universe. This will 
give a{p{t),x.) as a function of p{t) and cr, (x). Just before Co- starts to become significant 
at some epoch ^st.a.rt, this will give some relation CT(/9(tstart), x) = (7((T,(x)). We can then 
replace Eq. (pS)) bv [Hsj 

C.(x)~5A^(.g(a,(x)),0,). (42) 

The derivatives in Eq. (15(7)) are then evaluated using d/da^, = {dg / da^)d / dg. 

As long as a exchanges no energy with its surroundings, its evolution will be given by the 
field equation 

a + iHa + dV/da = 0. (43) 

This is a good approximation in almost all cases. (An exception within the original curvaton 
scenario will occur if a loses energy by preheating soon after it begins to oscillate [29|.) 
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IV. IMPOSSIBILITY OF CURVATON-TYPE PREHEATING WITH THE 

QUADRATIC POTENTIAL 

A. Inflation 



Now we consider preheating with the quadratic potential ([T]). We assume that the po- 
tential holds also during inflation, with the first term dominating so that V — to^0^/2. 
Under that assumption, we are going to show that curvaton-type preheating cannot occur, 
because the small value of g required to make x light during inflation is insufficient to cause 
significant preheating. In Appendix B we show that the same is true if the potential during 
inflation flattens out corresponding to what is called inflection-point inflation. 

With the potential y = m^0^/2, the inflaton contribution to ( gives V,^^ ~ m^/Mp. With 
= 10~^°M| this contribution dominates, which is allowed by present observations. We 
will take the value = 10^^'^ Alp as an order of magnitude estimate. (As will become 
clear, curvaton-type preheating is excluded even more strongly if m is lower.) Inflation ends 
when Mp making > M|. Therefore, x is hght only if g'^ < m^/A/p 10"^°. We 
write 35 = and require g^ ^ 1. 

During inflation the unperturbed fields evolve according to the slow roll approximation: 

3H^ = ~m^{(j))(t> (44) 

m2(0) = + g\^ = (^1 + g^^^ ~ (45) 

3Hx = - -gWx, (46) 

with 3MpH^ ~ m^0^/2. We assume x ^ Mp. The evolution of (f> is hardly affected by x, 
and well before the end of infiation N(t) — — ^^gr-, where N(t) is the remaining number of 

e-folds of inflation. Using this expression one finds xi^) [N{t)]^S'i, which means that x is 
practically constant during infiation. 

Inflation ends soon after the slow-roll approximation fails, which in turn happens when 
H{t) ~ TO. After inflation (j>'^ oscillates, and xi^) =oc exp[—g'^N{t)], where now N{t) is the 
number of e-folds after the end of inflation. This behavior persists until cf) decays, or the 
mass of X becomes signiflcant. These events must occur in time to create radiation with 
T ~ MeV, which means roughly N{t) ^ 10^. Assuming g^ < 10~^ the evolution of x is 
again negligible. We can therefore set x{t) = X*- 

At this point we mention the issue of radiative contributions to the potential, generated by 
the coupling g^ . According to the Coleman- Weinberg formula, the one- loop ^ contribution 
is 

Moop=g':^ln(0/Q), (47) 

where Q is the renormalization scale. To minimize higher order contributions one should 
choose Q ^ (j)^, and then the requirement that t^oop and V["^p are negligible corresponds to 
5^^ Aixm/cj}^ ^ 10"^. With supersymmetry the loop correction is reduced and it looks 
as though the small g that we are invoking is safe. Unfortunately though, this estimate of 
the loop correction is not self-consistent because the renormalization scale must be below 
the ultra-violet cutoff of the effective field theory, which in turn must be below Mp since we 
are assuming Einstein gravity. To bring the loop correction under control one would have 
to embed the model within a larger theory. We proceed on the assumption that the loop 
correction is negligible. 



B. Creation of x particles 

In the literature, the potential ([1]) is usually considered with not too many orders of 
magnitude below 1. Then the oscillation of (j) after infiation leads to the copious production 
of X particles. This occurs through what is called parametric resonance, and it rapidly 
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drains away most of the oscillation energy corresponding to what is called preheating. We 
are going to show that is too small for parametric resonance. Then we will show that 
the variation of around the end of inflation still creates some x particles, and go on to 
estimate the contribution of this creation to Cx ■ 

It is sufficient to consider narrow resonance, which requires the condition q = g^cjp' /Am? <C 
1. (The alternative regime of broad resonance, g ^ 1 requires a bigger value oi .) Narrow 
resonance occurs if there is time for (j) to undergo some oscillations, while the wavenumber 
k/a(t) passes through the narrow band 

l-q{t)/2<-^<l + q{t)/2. (48) 
m a(t) 

The time for the passage to occur is At ~ q/2H , so narrow resonance requires mq/H ^ 1. 
When the oscillation begins at iJ = ~ m, this corresponds to g^ ^ 1. Afterwards At 
decreases like a~^^^. It follows that narrow resonance requires g| ^ 1, which contradicts 
the light field requirement g| ^ 1. Hence the light field requirement is incompatible with 
parametric resonance. The estimates of non-gaussianity from preheating in [l9l . [20| are 
therefore inconsinstent, because they invoke both parameteric resonance and the assumption 
that X is light. 

The preheating process operates on scales that never leave the horizon. It occurs because 
the evolution of the mode function Xk is non-adiabatic. At some level, the evolution of the 
mode function will be non-adiabatic even if g^ is too small for preheating, which means that 
X particles will still be created. Let us see if this creation can give a significant contribution 
to C- 

We work with fk = axk and conformal time 77. Denoting d/dr/ by a prime, it satisfies [2| 
+ [k' - X^irj)] h ^0, X^^ a" /a - g^^a^ (49) 



To calculate the occupation number one imposes the early time condition \/2kfk = 
exp{—ikr]) . The late time behavior is then 

V2kfk = ake-''''^ +(3ke''\ (50) 

with |afep — — 1, and the occupation number is Uk — lA-P- 

We set ry = at the end of inflation. Well before the end of inflation, 77 = — l/aiJond 
and well afterwards rj = 2/aH. In both of these regimes, jryj ^ l/flcnd^^cnd where the 
star denotes the end of inflation, and we have X'^ ~ a" /a ~ 2/77^. We conclude that X^ 
has a peak centered on = with width of order l/ctcnd^^ond, arid that within the peak, 
X'^ ~ a" /a ^ l/(acnd^ond)^- As an approximation we may write 

a"/a- 2/(772 + (l/a,„difc„d)'], (51) 

which becomes accurate in the regime \r]\ ^ l/(acnd^^cnd)2- 

In the regime k » Oend-f^end, the square bracket in Eq. (|49l) is close to fc^ at all times, giving 
rife <C 1. In the opposite regime k <C Oend-ffend, which we are not considering at the moment, 
the vacuum fluctuation is promoted to a classical perturbation corresponding to occupation 
number nk ^ I. It follows that in the regime k > Oend^^cnd under consideration, particle 
creation occurs mainly at the bottom end of the range corresponding to fc ~ Ocnd-ffcnd and 
that nfc ^ 1 in this range. Assuming that these particles dominate we have 

3Px -Px--, r dkkW - i/eid (— ) ' . (52) 

(We present in an Appendix estimates of rife at fc ^ flcnd-ffond, showing that the contribution 
from this regime is indeed negligible.) These fall like a~'^ since we are taking x to be massless 
and the creation occurs at a ~ Ocnd- The total energy density is dominated by 

^ MlHl^^{a,^a/af. (53) 
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C. Estimating the contribution to 

The contribution (^^{x,t) is given by Eq. ([30]) with a — x- K is neghgible just when 
the smoothing scale leaves the horizon, which we are choosing to be a bit shorter than the 
shortest cosmological scale. This epoch is 50 or so e-folds before the end of inflation, and 
we denote it by a subscript 1. 

During inflation, (x, t) will increase slightly, because the presence of x means that there 
are really two slowly rolling fields. Because we are choosing <^ 1 and x Mp, this effect 
is negligible. We are interested in the possible increase of (x, t) shortly after inflation, due 
to the creation of the x particles. We will denoting this increase by C^. It is equal to 6N, 
where N is the number of e-folds taking place while the x particles are being created. 

Using Eq. ^ to first order in 5P, we have 

aCx dN dP/dx.da 

- - ' , (54) 



dx* dx* Ja^^a P + P « ' 

where a subscript 2 indicates an epoch just after the end of preheating. Without the tildes 
and with the integration going from the initial epoch to a sufficiently late time, this expres- 
sion gives the complete quantity dCx/dx*- It has not been given before. Note that x* is 
evaluated at the initial epoch during the inflation as discussed after Eq. (|22p . 

To evaluate dP/dx*, we use Eqs. (|53)l and (|52|) for p and P. They depend on x* because 
the epoch at which inflation ends depends on the effective inflaton mass to, which in turn 
depends slightly on x* according to the expression H^^^ ~ to^(x^). Inserting this expression 
into Eqs. ((53)) and ([52]), and differentiating at constant p, we find 



- -in? 



dim'') 9 "^"n a 
Using d{m''{x*))/dx* = g^X* and Eq. ([Si]) we find 



dN g'^x* _ „2'7^^X* 
c 

From these expressions we find 



dx. ' M2 - ^5 ■ ^^^^ 



P,/.,..„^,.(ZL)*(i.)%<(ZL)'.,o-» (57, 



Assuming x* ^ ^^cnd the contribution of (Sx*)'^ to the non-gaussianity parameter is given 
by Eq. §^ 



6 / \ 2 / X 6 

-30 



10-^^ (58) 



If instead x* -f^cnd the result for f^^ is even smaller, given by Eq. (|571) . 

Of course this is not the only contribution to /nl- The contribution N-^Sx* will contribute 



54j because the bispectrum of dx* is not exactly zero |45|. The dominant contribution 



though will be from the inflaton perturbation, giving [45|,|46| |/nl| 10 



V. MODULATED PREHEATING WITH THE QUADRATIC POTENTIAL 

To generate significant preheating with the quadratic potential ([T]), we need g' ^ 10^^°. 
This makes x heavy during infiation, and for preheating to generate a significant contribution 
to ^ we need to introduce a third field a. The contribution of cr to the potential should be 
such that (i) a is light during infiation and (ii) its value contributes to the effect value 
of g so that we have g{(y*). Then the perturbation Sa on cosmological scales will give a 
non-adiabatic pressure perturbation between the end of inflation and reheating, which will 
gen erate a contribution to This scenario is similar to the modulated reheating scenario 
[iSl - llTt . where the decay rate T{a*) depends on a light field x- For comparison, we recall 
the modulated reheating result, using the 6N formalism [stJ. 
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A. Modulated reheating 

For modulated reheating it is a good approximation to take the decay as instantaneous. 
Then vanishes before decay and is constant thereafter. It can easily be calculate because 
the dependence of the local scale factor is known both before and after decay [Tsl - flTl . |37| . 
Before decay there is matter domination with a cx The decay occurs when H = T, 

therefore 

(59) 

This gives cx T"" with a ^ 1/6. Then 



This gives [37 



dN 1 d^N 1 



= -aT'/r, = ~aV"lT + a [V /v f , (61) 



and 



a 



p/2 



(62) 



For illustration, one may take a simple possibility [17|: 



=.9^(1 + ^1, g{a^)^gll + ^j^\::.g. (63) 



with M some mass scale. Then 



/nl ~ (1-:^). (65) 



These expressions allow Co- to dominate and they allow /nl to be much bigger than 1. 



B. Modulated preheating 

Preheating generates a contribution to ( which changes continuously during preheating, 
and continues to change afterwards for as long as generates a significant non-adiabatic 
pressure perturbation. We have therefore 

Ca=SN = 6Np,, + ,5iVafter, (66) 

where Np,-c is the number of e-folds of preheating and A'aftcr is the subsequent number of 
e-folds up to an epoch when the non-adiabatic pressure has become negligible. 
Given N{a^,) we have 

«, ^ (67) 

„dN ,2d^N , , 

Inserting these formulas into Eqs. ([1]) and ([37)) gives the preheating contribution to the 
spectrum and non-gaussianity of the curvature perturbation. 
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FIG. 1: /nl in the plane of x/M and M in case of = lO"'^. We plotted /nl = -10"^ -0.1, -1 
and —10 from top to bottom. The boundary Vc„ < Vc.o\,s ( '--^ 3 x 10^^) is also plotted by the solid 
line. 



For simplicity, we will assume that p^^ has almost redshifted when decays, so that the 
non-adiabatic pressure has become negligible by that time. Since cx after preheating 
is over, this gives 

7Vaftor(x,i)- Jln(^^), (69) 

where the subscript 2 indicates the end of preheating and t can be taken as the time just 
before </> decays. We then have 

3^ = (70) 

og P4>2 og 

In general, an accurate calculation of and iVprc requires numerical simulation. Here we 
will just use analytic estimates which should be good in the regime 10"!° < g^ < 10"^. 
In this regime, Pv2 ^ P<p2, with the equality attained for ~ IQ-^. 

According to [5|, the x particles are non-relativistic during preheating. Hence p-^ as well 
as p0 is matter-dominated during this era and the total energy density is proportional to 
a~^. According to [s'], preheating ends at the epoch given by 

imt,r-(^X-(^y=gl (71) 



^cnd 



where <2 = gv^M-p /2>w? and g^ = ff/lO ^. Let us check that the x particles are indeed 
non-relativistic up to this epoch. Their effective mass-squared is 



ml{t) = 52^2 _ g^M'^{a,,,d/a{t)f , (72) 
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FIG. 2: Same as Fig.[I]but for = IQ-* 

where the star denotes the end of inflation. According to 5], the momentum of the x 
particles in our regime is Q k/a ^ ^^gMpm{acnd/<^)- It follows that the x particles are 
non-relativistic until 

(a/aond) ^ g— = .95, (73) 

m 

which is after the end of preheating. 
From Eq. ([7T|) we get 

dg 3g' dg^ ig^' ^ ' 

To handle A^aftcr we need The total energy density at the end of preheating 

follows from Eq. ((7T|) : 

Pond ^ M|m^(aond/a2)^ ~ / g"^ ■ (75) 

Section IXA of @ nicely gives the expression of n^{t) for t < tend- Multiplying it by m^{t) 
to get , we find at the end of preheating 

pM = Ag-^e^\ (76) 

with 

35/2 X lO-^m^ 

A = — ^ — : (77) 

B = 2x (87r)i/2 X O.UMp/m ~ Mpg/m. (78) 

The difference between Eqs. ([76]) and ((75|) gives p'^'^ig), from which we can calculate 
9iVaftoi/9g and S^TVaftor/^g^ as follow, 

dg 3(p™d„^c„d) [ g Px )^ ('^) 
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dg 



after 1 

^ 3 



-1 



(^ond _ pOnd)2 



2p- 



-Bpl 



p2 _cnd 



(80) 



Adding these into dTj), we obtain the total dN /dg and d'^N/dg'^ 



Except near the top of our range 10 ^'^ < g^ < 10 the ratio = p'^'^/pend is much 



less than 1. To first order in i?^ we find simply 

ON 

— ^ -(Mp/m)i?, 



-WRy 



IQi^i?, 



(81) 
(82) 



Taking g{x*) to be given by Eq. ([63|) we can calculate T'^^ and /nl from Eqs. ([32|) -([37 |) . 
In the regime where either of them may be significant, the loop contributions dominate. 
Setting In(fcL) = 1 we plot the result in Figures [1] and [H (The exact result was used instead 
of Eqs. (|8T|) and (|82)). which is needed for Figure [TJ) We see that the observational bound 

/i 



NL 



> 



-10 can be satisfied only if Vq^/Vq is very small. 



VI. MASSLESS PREHEATING 



In this section we consider curvaton-type preheating with the potential ([2]), usually called 
massless preheating. Curvaton-type preheating, corresponding to a contribution is 
indeed possible for this case, because there is a region of parameter space which makes x 
light during inflation. Also, the inflaton trajectory can be in practically the direction. 
Then we we deal with single-field inflation, and the contribution ^-^ becomes significant only 
after inflation is over, ie. the preheating contribution to is the dominant one. 

With the potential V = A(/>'*/4, the inflaton contribution to C gives [37|, |55| Vi^^ = lO^A. 
To agree with observation it should not dominate [s^ which requires A < 10~^^. To obtain 
significant preheating one requires roughly ~ A. These small couplings beg explanation, 
and as with the quadratic potential there is also the issue of keeping radiative contribution 
to the potential under control (though the naive estimate 14oop ~ ff*'/'^ ln('/>/Q) makes that 
contribution tiny). We proceed without addressing these issues. 



A. Inflation 



To keep things simple we will pretend that all cosmological scales leave the horizon 55 
e-folds before the end of inflation, and denote their values then by a subscript 55. Following 
[23] we focus on the case = 2 A (corresponding to a supersymmetry) . 

Single field inflation occurs if dV/dcf) ^ dV/dx, which corresponds to x ^ 0- We assume 
this condition after the biggest cosmological scale leaves the horizon, which occurs when 
(p ~ lOMp. Then, during inflation and for some time afterwards, 

4> + 3iJ0 + A03 ^ 0, (83) 

with = A</>'^/12Mp during inflation. Inflation ends at the epoch 0ond — Mp. During 
inflation, t]^^ = 4Mp , so that x is light except near the end of inflation. 
In the slow-roll approximation, the evolution of the fields is given by 

^ iMl (84) 



H 3iJ2 

X _ g^4''^x _ ^g^M^x 



(85) 



H 3iJ2 3 A (/-2 • 
Well before the end of inflation, these give 

</>(i) ~ v/8iV(I)0end X{t) = (8iV(t))3'/6^Xcnd, (86) 
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where N{t) is the remaining number of e-folds of inflation. The earhest epoch of interest is 
the one when the scale k — Hq leaves the horizon, corresponding to N{t) ~ 55. Then 

055 ~ 21(/)cnd X55 - SXend- (87) 

The corresponding Hubble parameter is i?55 ~ 440i/ond- An accurate calculation [s^l, taking 
account of the failure of slow roll near the end of inflation gives X55 ~ lO^Xcnd- 



B. Preheating in an unperturbed universe 

The field equations are 

4> + 3H(j) + X(j)^ +g'^(l3x'^ = fl-^V^f/) (88) 
X + 3Hx + g^^\ = a-^V\. (89) 

In the regime x* ^ 'P*! the first order perturbation of Eq. (IMl) gives 

Xk + 3i7xk + g'0\k + ik/a)\k = 0. (90) 

In the early stage of preheating, <j) oscillates according to Eq. (|83| . This corresponds to a 
time-averaged P = p/3 (not to P ^ p as for the quadratic potential). The broad resonance 
band of k is time-independent (not redshifted as for quadratic inflation) and depends only 
on the combination g^/X. For cosmological scales to be in the band, it has to extend down 
to k = 0. That is the case for 1 < g^/X < 3, which is satisfied by our choice g^ — 2X. There 
is no significant preheating for g^ < X. 

As we noted earlier, C ~ is forbidden in this model. In the opposite case ~ the 
spectral index is given by Eq. (|4ip . In the slow- roll approximation (|86p with x — this 
gives'^"* 

This too is forbidden by observation with our chosen value g^ ~ 2X but it might agree with 
the minimum value g^ = X. Otherwise we would need to assume both Q ^ C (-^ ^ 
requiring a third light field to give the dominant contribution to (f>. 

To arrive at an estimate of C^, we will invoke (27j which computes the number 
A^pre(Xend, 0cnd) of e-folds of preheating in a universe that is unperturbed at the end of 
inflation. To be more precise, iVpi.e(xcnd, 0cnd) is the number of e-folds from the end of 
inflation (when H = —H^) to an epoch of fixed p soon after the end of preheating which is 
around five e-folds later. 

The calculation is done within a comoving box, whose size at the end of inflation is 
20/ffend, with 7Jend/27r = 10~^ X -y/87rA/p.^^ The lowest wavenumber is then k = 
(7r/10)(a_ff)cnd, which means that all k come inside the horizon soon after preheating begins. 

Two methods of calculation were used, which give similar results One of them 

uses two stages [s^l- The first stage, lasting for about two e-folds, is done in Fourier space. 



We fixed the coefliicient 0.02 by taking slow roll to be valid right up to the end of inflation, with 4>'^^^^ = 8Mp 
which corresponds to ending inflation at e = 1. A numerical calculation handling the departure from slow 
roll at the end of inflation alters this estiamte by only a few percent [57l . 

At the final epoch, the cosmic fiuid has equation of state P/p = 1/3, except for the contribution of 
the X oscillation. The latter is small, so that p(t) is smooth except for a small oscillation. This small 
oscillation is averaged over so as to avoid the small oscillation in A^pro (Xend 1 </'ond ) a function of Xendi 
that would otherwise be present. In order not to affect the final outcome, represented by the smoothed 
function Nii{xend><Pcnd defined below, that oscillation should have period <^ fl^cndi but we have not 
checked whether that is the case. 

Note that their Mp is %/8n times ours. The value of -ff^nd is about the one that would make (^^ equal 
to the observed value. In reality, observation requires and therefore -ffend to be somewhat smaller but 
that would hardly the following analysis. 
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The k = mode is evolved using the non-linear equations (|88p and (|89|) (with the right hand 
sides zero). The k ^0 modes are evolved analytically using Eqs. and Before each 
such mode leaves the horizon, it is treated quantum mechanically so that Eq. (f83| applies 
to the mode function. In this way, the vacuum fluctuation present before horizon exit 
becomes a classical perturbation after horizon exit."*^^ Up to some /cmax there is parameteric 
resonance generating classical perturbations from the vacuum fluctuation. The second stage 
is a numerical simulation of Eqs. ([55]) and (|89l) using a new code Hlattice. The other 
method does the numerical simulation from the beginning, using an update of the existing 
code DEFROST [581]. 

The result using DEFROST is presented in Figure 1. of [13], with ^end fixed at the 
value corresponding to e = 1 and with Xend in the range 10~^Mp <yend ^ 10~**Mp.#^ 
In accordance with earlier (but not sufficiently accurate) work [y, [2y|, it is found that 
-^pre(Xend, (/"end) is not smooth function of Xond, but rather exhibits sharp spikes. There are 
big spikes having a spacing of order i?end with smaller spikes in between them. The width 
of the spikes (not visible on the plot) is of order 10~^i/ond fs?]. 

As explained in 27], the spikes occur because the evolution of the zero mode {0(t), x{t)} 
is given, until a fairly late stage, by the unperturbed equations (Eqs. (ISS)) and (ISQ)) with the 
right hands sides zero). Each spike corresponds to a value of {^icnd, Xcnd} that makes x{t) 
anomalously small at some epoch . 



C. An estimate of (^^ 

We are trying to calculate the contribution C^, of Sx* to C. It is defined by Eq. (^5)) and 
we would like to arrive at the approximation Eq. (15(11) . In [23|, it is proposed that this is 
done by following procedure. First smooth Np-^dxcnd, '/'end) (as a function of Xond with fixed 
9^end) by a gaussian window function with variance 

al = (i/ond/2^)' ln[i?(ai7)e„d], (92) 

to obtain a function Nj^^Xend, 4'end)- (As the notation indicates, cr^ is the mean square of 
f^Xond, within a region of size R.) Then write 

Cx(x) = 5Affl(Xc„d(x), (/-end). (93) 

The smoothed function iVfl(xcnd, 0cnd) is shown in Figure 1 of [27| . 

We will consider the validity of this procedure in the next subsection. Adopting it for 
the moment, we can use information given in [27| to estimate V^^ and /nl- In the range 
X* ^ 10~^Mp, the smoothed function Nn is quadratic: 

Nr = Uxl /x = 2 X IOVSttMI. (94) 

Let us suppose first that x* (the spatial average within the minimal box) vanishes. (As 
we noted in Section fill Dl that is unlikely but it cannot be ruled out.) Then Eqs. and 
give 

Since i?*/27r < 10^'' Mp we see that x* is giving a negligible contribution. 



Modes with k <C (a-ff)cnd are classical already at the end of inflation, but the small box size practically 
excludes such modes, justifying our statement that the simulation starts with an unperturbed universe. 
An arbitrary constant is subtracted from A'^prc(Xend, '/'end) to make it of order 10~^, and the result is 
somewhat confusingly labelled as SN. 
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This calculation of /nl may be compared with the one reported in [25|, which also took 
the unperturbed value of x* to vanish. For they invoke the usual formula of first order 
cosmological perturbation theory, which as is well known is manifestly the same as the 
6N formula ([25)1 . But for Cx they invoke second-order cosmological perturbation theory, 
stopping the calculation at an epoch when just a few oscillations of (f>{t) have taken place, 
so that parametric resonance still applies. This should in principle give a valid result for 
C,x at that epoch, though one might be concerned about the complexity of the second order 
equations'*^ The final step would then be to calculate the spectrum and bispectruni of 1^,^,, 
to obtain 7^^^ and /nl- That step was not however performed in [25j . Instead, a formula 
/nl ~ C2/C1 was invoked with ^2 = Cx Ci = C^- But this formula applies only if 
is of the local form, ie. if C2(x) = (3/5)/nlCi (^)- In. the present case where Ci and (2 are 
uncorrelated it obviously cannot apply. Neither does it apply if the quantities are replaced 
by typical magnitudes so that it becomes /nl ~ T^l^"^ /'Pc^,- The correct estimate [H^ is in 
fact /nl ~ ^Cx^/^C*' corresponding to Eq. (1571) . 

So far we took x* to vanish. Now assume instead a more likely value x* ^ -f^*- Then 
Eqs. ((33)) and ([36]) apply, giving 

/nl 

According to Fig. 1 of [l^l, the quadratic approximation holds only out to x* ~ 10~^Mp, 
which means that Vq^ and /nl are again negligible. We conclude that if the unperturbed 
X* is small enough for Eq. ([M)) to apply, the curvature perturbation in this scenario is still 
negligible at the end of preheating. 

In the range 10"^ ^ x*/Mp < 10^'', the presented in Figure 1 of [l^l is an oscil- 
lating function of Xcnd- As one would expect, the oscillation is slow so that Nr is still a 
smooth function of Xcnd over an interval A^ond ^ ^cnd- Therefore, the usual approximation 
Eq. ()54p should still be good. The corresponding spectrum and bispectrum have not yet 
been calculated, and the calculation of N{x*) has not yet been extended to larger x*, but 
the smoothness of Nn and hence the validity of Eq. (IMl) should still hold. Observational 
consequences of its failure are mentioned in in [13] . 



4 X 10" 



= 5 X 10" 



lO^X. 
Mp 

10^ 



2n 



27r 



(97) 
(98) 



D. Justifying the estimate of 1^^ 

To formulate an approximation within which Eq. (1931) can be justified we proceed in 
two steps. First we formulate an approximation within which A^pre(Xend(x), ^end) can be 
regarded as the number of e-folds of preheating at a given location. Then we justify the 
smoothing procedure that gives iVi?(xond(x), 0cnd)- 

The first step requires the existence of a scale Lhom, which is big enough to be well 
outside the horizon at the end of preheating, yet small enough that the variation of 
A^pro(Xcnd(x), 0cnd) IS negligible over a distance Lhom- Let us consider these criteria in 
turn. 

Taking a oc t^^'^ during preheating (corresponding to roughly p ~ P/3) and assuming 
A^pre ~ 5 e-folds of preheating, the first requirement becomes 

ihom » e^^-H-^\. (99) 



Analogous equations presented in a two-field inflation model lead to a result |36l| that disagrees with that 
obtained from the 5N formalism [35j . As the 5N calculation in that case is very straightforward, the 
discrepancy in that case presumably indicates an error in the perturbative calculation, caused by the 
complexity of the perturbative equations. 
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For the second requirement, we can take 'negligible' to mean 'much less than the height 
^-^peak of a typical peak in -/Vpi.c(xcnd)- The variation of -/Vpi.c(Xcnd(x)) within a distance 

Lhom is 

A7V^^^|VXo„d|ihom. (100) 

«Xond 

The maximum of dA^pro/dXcnd is achieved near a peak, and is of order AA^poak/Axpcak where 
Axpcak is the width of the peak. The requirement AA^ ^ AA^pcak is therefore 

|VXend|iliom < AXpoak- (101) 

Using Eq. (jl3p to estimate | Vxcnd | , this requires Xcnd to be smooth on a scale R satisfying 

R » HcndLhom/ AXpoak- (102) 

Using Eq. ([M)) this becomes 

R{aHUd » _^e^-. (103) 

^Xpcak 

Using the estimate Axpcak lO^^ifcnd, with A^pro — 5, the right hand side is of order 
which means that Xcnd should be smooth a scale of order (aH)^^^. To achieve this, 
we have to drop the perturbations of Xcnd that are generated during the last 10 or so e-folds 
of inflation. There seems to be no justification for this, but let us anyway move on to the 
second step. 

We now have an approximation which makes iVpi.c(x) the number of e-folds of preheating 
at a given location. But Np,-c{x) is s not yet the quantity whose perturbation is (^^. The 
reason is that the latter is, by definition, smooth on the chosen scale R whereas iVpi.e(x) is 
not because of the narrow spikes. Before using the formula = SNp,:^, we must smooth 
iVpre(x) on the scale R. We are now going to argue that the result of this smoothing will be 
approximately 7Vj^(x). 

Smoothing on the scale i? at a given location x means that the quantity is replaced by its 
average within a sphere of radius R. At a random location within such a sphere, Xond has 
a gaussian probability distribution with mean square cr^ given by Eq. (j93p . If Xcnd(x) had 
negligible correlation length (one much less than R) the distribution at different locations 
would be independent and we would immediately obtain the desired result. Unfortunately, 
the correlation length is practically infinite because the spectrum of Xcnd(x) is practically 
flat at long wavelengths. We therefore have to proceed differently. 

We need some notation. For any function 5 of a variable y, let us write the smoothed 
quantity as 



g{y,ay)= / dy'g{y')W{y' - y,ay), (104) 

-/ —00 

where the wmdow function satisfies"^ 

/oo 
dy'W{y' -y,ay) = 1 (105) 
-00 

W ~ i\y' - y\ » ay) (106) 
W ~ const {\y' ~y\<^ay). (107) 

The window function is usually taken to be either a theta function (top hat) or a gaussian. 
In any case, ay is usually taken to be a constant. Then the convolution theorem shows that 
smoothing kills the Fourier components of g for wavenumbers k ^ ay while hardly affecting 



Equivalently, one can replace the third condition by ly ~ 1, and divide Eq. ||104|I by the left hand side 
of the first condition. 



21 



them for k ^ ay. The same wih be true even if <jy depends on y and/or y' provided that 
its variation is neghgible in the regime \y' — y \ <^ Oy. 

Let us denote iV(xcnd, '?^'cnd) simply by N{x), and iV_R(xcnd, 0cnd) by Nb.{x). We have 

Ndx) = J dx'N{x')W{x' - X, <^x), (108) 

with given by Eq. ([M)) . We want to show that iVi^(x(x)) is smooth on the scale R. Let 
us first assume that Vx is practically constant within a region of size R. Then, choosing 
the z axis in the direction of Vx, 

X(x')-x(x)^^(z'-z), (109) 

where x' has coordinates {x,y,z'). We then have 

7Vk(x(x)) ~ J dz'N{x{^'))WRiz' - z, a,), (110) 

where 

WR{z'-z,a,) = W{z' -z,a,){dx{-y^)/dz)-^ (111) 

fT,(x) ^ a^{dx{^)ldz)-\ (112) 

From Eq. (jl3p we learn that crz(x) ~ i? at a typical position. Therefore, as required, 
Affl(x(x)) is just the original function A^(x(x)) smoothed on the scale R. 

We have assumed that Vx is practically constant within a typical region of size R. To 
see whether this is reasonable, we can estimate the mean squares of the second partial 
derivatives didiX, whose Fourier components are —kikj times those of x- Analogously with 
the calculation leading to Eq. (IT5)) . we find that the typical fractional variation of Vx within 
the region is of order 1. To handle this order 1 variation, we could use curvilinear coordinates 
with z constant on each spatial slice of constant Vx, and our conclusion about the Fourier 
components of N]i{x{x)) would still hold. We have not considered the case where Vx has 
very large variation, since it will apply only to rare locations that should not affect our 
conclusion. 

The main source of error in this prescription for (x) is the neglect of the contribution 
to (5xond, that is generated from the vacuum fluctuation after R leaves the horizon. We have 
no idea how to estimate this error, except by performing a numerical simulation in a box 
with size rather bigger than e^(27r/(a_ff )cnd) that would support those modes. Even with 
such a simulation in place, there is at least in principle the fundamental problem, that the 
result for Cx(x) would depend on the chosen realization of the small-scale perturbations. 
One may hope that this dependence is small. If it is not, the derivative would become a 
stochastic quantity, and instead of T'^^ = ^x^x would have V]^^'' — {N^)V^. Similarly, 
in calculating say one would have to replace N^^/N"^ by its expectation value. These 
expectations values would have to be calculated using the numerical simulation. 



VII. CONCLUSION 



In this paper we have carefully addressed some issues, that arise when one uses the 5N 
formalism is used for preheating. Then we saw how things worked out in a couple of specific 
cases. 

The investigation seems worth continuing, in a number of directions. Several of the 
the curvaton-type preheating scenarios listed in the Introduction have been explored only 
with cosmological perturbation theory. It seems worthwhile to look at them also with the 
5N formalism, especially in cases where second order perturbation theory has been used. 
The other paradigm, modulated preheating, would also be worth exploring, starting with a 
numerical simulation for the simplest setup using the quadratic potential of Eq. ([T|). 

All of this assumes that the curvature perturbation is generated exclusively by one or 
more scalar field perturbations. As has recently been realised [59l - l62| . it might instead be 
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generated, wholly or partly, by one or more vector field perturbations. A smoking gun for 
such a setup would be statistical anisotropy [60l - l62j . The vector field possibility has so far 
been explored only with the generation of ^ taking place through the curvaton mechanism 
before a second reheating during inflation [6l| or at the end of inflation [g^]- It 

is clear that one could implement modulated preheating using a vector field perturbation, 
allowing one or more of the parameters to depend on the vector field in the spirit of (60l |. 
One might also implement curvaton-type preheating using such a field, provided that the 
preheating can create a vector as opposed to a scalar field. It is not known at present 
whether that is the case. 
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Appendix A: Occupation number of x quanta 



In this Appendix, we discuss the solution of Eq. (H^ , 

f'l^k^h^Xh, X = -g^<i?c? + a!' la, (Al) 

and estimate the occupation number of the x quanta. We work in the regime fc 3> a*i?* 
where the star denotes the end of inflation. According to the estimates after Eq. ([50|) . AT/fc 
is a small perturbation in this regime, which we treat to first order. We write 

h^h.^ + ^h. (A2) 

where /fe_o — exp(—ikr]) / ^/2k is the unperturbed quantity. When this equation is substituted 
into Eq. we get a following equation at the first order of the perturbation, 

Sf'k + kHfk = Xk{ri), Xk = Xfk,o (A3) 

To solve Eq. (jA3|) , we use the Fourier transformation method which requests variables to 
be expanded into their Fourier modes, 

Sfkiv) r di. 6Uu) e'-^ (A4) 

and 

Xkiv) r dLo Xk{Lo) e^"". (A5) 



27r „ 

In reverse, we can define the inverse-Fourier transformation of by 



XfeM= / dT^' Xkir,') e-'^^' . (A6) 



Substituting Eqs. (|A4p and (|A5p into Eq. (|A3p . we immediately get a relation between the 
Fourier modes. 

Therefore we finally have 
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The integrand has two poles, and to impose the condition 6fk = at — — oo we deform 
the path — oo < w < oo so that the poles lie below the path as shown in Fig. [3] (a). At 
7/ — > — oo we obtain the required initial condition by closing the contour in the lower half 
plane. To obtain instead the behavior at 77 — ^ +00 we close it instead in the upper half plane 
shown in Fig. [3] (b) which gives the behavior Eq. (15(11) with 

= pk) /I dve-^'^'^iv) (-ffV^a^ + ^) . (A9) 
Using the approximation (j51l) we find a standard Fourier transform, which gives 

\/3k\^ =n^{a,H^/kfe-^\''\ {k:^a,H,). (AlO) 

Putting this into Eq. (j52[) we find that it gives a negligible contribution to , as advertised 
in the text. One might be concerned that the exact result for |/3fcp might fall off more slowly 
at large k giving a significant or divergent contribution to |/3fep, but the following argument 
should allay such concern. Since X is infinitely differentiable, integration by parts shows 
that 

Pk = / drye-^^'^"^, (All) 

-2fc(2fc)" ' dry" ' ^ ' 

for all n > 0. For low n, it is reasonable to extend the argument leading to Eq. ([5T|) . to 
arrive at an estimate 

v(") - . ^ \ ( ^^'>\ 

This gives 

|/3fe|2 < (a,F,/fc)2+", {k:^a,H,). (A13) 

Using this as a reasonable approximation for n = 3, we again find a negligible contribution 
to px- 



Appendix B: Preheating after inflection point inflation 

In the main text we have studied the preheating only in chaotic inflation models, where 
the potential V = \m?(j)^ generating the oscillation holds also during inflation. In that case 
(the value at the beginning of the oscillation, or order its value at the end of inflation) 
and m are given by 0* ^ A/p and m? < 10^^*^, with the latter inequality saturated if C,tp is 
to be a significant fraction of the total. Also, the Hubble parameter at the beginning of the 
oscillation, given by 



(Bl) 
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is of order m. 

In a different class of models, known as inflection point models, the potential flattens out 
at > so that its first and second derivatives are close to zero during inflation. Then 
0* and m are independent parameters and ff* given by Eq. (|B1[) is much smaller than m. 
We now analyze preheating in these models, assuming the interaction ^g'^cjy^x'^ ■ Of course 
the viability of that interaction and the possible identity of x should be checked within a 
particular setup. 

Here we see what happens if is smaller. The discussion may be relevant for of inflation, 
either in the context of MSSM inflation IgJi which have 0o ~ TeV or in the context of 



colliding brane inflation [6J] where (/)o might be anywhere in the range TeV ^ ipo ^ Mp. Of 
course, one has to check within a specific model if the interaction term g'^cfy^X^ is consistent 
and reasonable. 

The Hubble constant at the end of infiation is 

When we consider the preheating which would be induced by an interaction term such 
as the second term in Eq. ([1]) the effective mass of x is given by — g'^cj}?). We can 
parameterize the amplitude of the oscillation term in the Mathieu equation as [3| 

9 = 5^4- (B3) 

First of all, we shall consider a condition for lightness of x during infiation, <C H^. 
Thus, it is found that 

.««^, (B4) 



which means 

lQ-3of^y_ (B5) 



Here only in case of MSSM inflation, ^ TeV. 

Next we discuss the condition for successful parametric resonance. As was discussed in 
Section llV Bl the parametric resonance occurs if k/a can stay in the resonance band for 
a sufficiently- long time to oscillate many times. The time interval when k/a passes the 
resonance band is given by At ^ q/H^. Thus that condition is represented by z = m^Ai ^ 
1, which gives 



=5 — ^ > 1- (B6) 
11^ nil 



2 

mi 



Then we have 

(poMp 

This lower bound on g^ is much bigger than the upper bound required by the lightness given 
in Eq. (|B5I) . when (po <^ Mp as is the case for the A-term inflation. 
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